June 16, 2005 



Sphalerons in the NMSSM 

Koichi Funakubo a l , Akira Kakuto 6 ' 2 Shuichiro Tao a ' 3 and Fumihiko Toyoda 6 ' 4 

^ | ^Department of Physics, Saga University, Saga 84-0-8502 Japan 

£ 
3 



X 



^School of Humanity- Oriented Science and Engineering, Kinki University, Iizuka 

820-8555 Japan 



Abstract 



We study sphaleron solutions in the next-to-minimal supersymmetric standard 
\q \ model. We find that the boundary condition on the singlet field at the origin of 

the radial coordinate is of Neumann type, while the other boundary conditions 



are of Dirichlet type. The sphaleron energy takes almost the same value as in 



O the MSSM for wide range of parameters, in spite of the negative contribution 

^ \ from the cubic term in the Higgs potential. 

cl 

Ph! 1 Introduction 

Since the discovery of the sphaleron solution[l], unsuppressed transition of the baryon 
number at high temperatures has been recognized to play an important role in early 



universe. It yielded new possibilities of matter generation, baryogenesis through lepto- 
genesis and electroweak baryogenesis [2]. The latter requires that the sphaleron process 
decouples just after the electroweak phase transition, in order to protect the generated 
baryon asymmetry from washout by sphaleron process in equilibrium. This sphaleron 
decoupling condition is expressed as r sph (T) ~ Te~ E ^ JT < H(T), where E sph is the 
static energy of the sphaleron solution and H(T) is the Hubble parameter at temperature 
T. In the standard model with one Higgs doublet, this condition can be cast into the 
form of vc/Tc > 1, where Vq is the expectation value of the Higgs field at the transition 
temperature Tc, where the gauge symmetry is spontaneously broken. Thus, the knowl- 
edge of energy of the sphaleron solution in a model viable for electroweak baryogenesis 
is indispensable to qualify the model by examining the sphaleron decoupling condition. 
The sphaleron solution was originally found in the 4-dimensional SU (2) gauge theory 
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with one Higgs doublet [3]. The minimal standard model, however, cannot afford to gen- 
erate the baryon asymmetry at the electroweak phase transition with an acceptable Higgs 
mass. It was also pointed out that the CP violation in the KM matrix is insufficient to 
provide the chiral charge as the source of baryon number. Hence, one must extend the 
standard model for successful electroweak baryogenesis. In particular, models with more 
bosonic fields are expected to admit the strongly first-order phase transition. Among 
those models are the MSSM with a light stop and the two-Higgs-doublet model. The 
sphaleron solutions have been found in the two-Higgs-doublet model[4] and in the MSSM 
with finite-temperature corrections [5]. The model with one Higgs doublet and one singlet 
is among those extensions of the standard model for which a sphaleron solution has been 
found [6]. 

We recently found that the Next-to-MSSM (NMSSM), which contains a singlet super- 
field in addition to the MSSM, can have strongly first-order phase transition, with heavy 
stops and Higgs bosons whose masses and couplings are consistent with results of collider 
experiments so far [7]. Although the NMSSM has more parameters than the MSSM, there 
are more constraints on these parameters which are absent in the MSSM. For example, 
the electroweak vacuum is not always the absolute minimum of the Higgs potential. This 
is because the cubic terms containing the singlet scalar can easily generate minimum of 
the potential other than the electroweak vacuum. We observed that the vacuum condi- 
tion, which ensures that the electroweak vacuum be the absolute minimum of the effective 
potential, together with the spectrum condition on the Higgs scalars, restricts allowed pa- 
rameters in the model[8]. In order to determine the parameter sets which are suited for 
the baryogenesis, one must know the sphaleron solution and its energy in this model. In 
the MSSM-limit, where the vacuum expectation value of the singlet goes to infinity, the 
sphaleron solutions are expected to have the same profile as the MSSM. When the expec- 
tation value is of order of the weak scale, we expect new features of the phase transitions, 
but even the existence of the sphaleron solution is not obvious. Our aim is to find the 
sphaleron solution and to study its energy in the case of weak-scale expectation value of 
the singlet field. 

Since the sphaleron solution is a saddle-point configuration, it is difficult to find it 
by solving the full equations of motion. Instead, we usually adopt the ansatz including 
the parameter along the noncontractible loop in the configuration space [3]. The ansatz is 
constructed in such a way that the configuration corresponds to the vacuum at the loop 
parameter fx = and tt, while the configuration at fx = it/2 is expected to be the highest- 
energy one corresponding to the sphaleron. In order to ensure that the configuration is 
the sphaleron, one must check that the configuration does satisfy the full equation of 
motion and that the fluctuation spectrum around the configuration does contain only 
one negative mode. In practice, these procedure are complicated and has been done only 
in the one-doublet model [9]. Here we extend the ansatz including the noncontractible 
loop to the NMSSM with [/(l)-gauge sector being turned off, and write down the static 
spherically symmetric equations of motion for the configuration at fx = ir/2. Then we 
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solve the equations of motion and evaluate the energy along the noncontractible loop 
around the solution to ensure that it is the highest-energy configuration along the loop. 

This paper is organized as follows. In §2, we extend the ansatz of Klinkhammer and 
Manton for the configurations along the noncontractible loop to the NMSSM. Then we de- 
rive the equations of motion for the sphaleron configurations. Several numerical solutions 
for the equations are shown in §3 for parameter sets corresponding to the light-Higgs and 
the heavy-Higgs scenarios. Section 4 is dedicated to concluding remarks. The derivation 
of the asymptotic solutions, together with the boundary conditions, are summarized in 
the Appendix. 



2 Ansatz and equations of motion 
2.1 noncontractible loop 

As discussed in [1], any static configuration of finite-energy doublet scalar $(a;) defines a 
map from S 2 spanned by the spatial coordinates (9, 4>) to S* 3 = SU (2) which characterizes 

the field at r = oo as a unitary transformation of the vacuum (^j- A one-parameter 

family of such configurations connecting vacua, which cannot be contracted to a point, 
is constructed by realizing the family of maps as a map from S* 3 to S* 3 . Among such 
one-parameter families, that of the highest symmetry is expected to form the least energy 
set of configurations. The configuration with the highest energy along the parameter will 
be a saddle-point configuration with one negative mode. 

The Manton's ansatz for the noncontractible loop in the one-doublet model is 

Sfor.M) = ^{(l-%))( e _4o S/U ) +/l(r)f/( ' U ^' 0) (l)}' (2 ' 1} 

A(/.,r,M) = --/W^^lt/- 1 ^,^), (2.2) 
9 

where 

jjt q a\ _ /e* M (cosyU — isin/icos^) e 1 ^ sin /i sin 9 \ (2 3) 

^' ' V — e -1 ^ sin // sin e~ v (cos/i + isin,ucos#) / ^ ' ' 

Here /i e [0, n] and the configurations are reduced to the vacuum at /i — and ir. As for 
the NMSSM, we adopt the following ansatz for the doublet Higgs fields, $d and $ n , and 
the singlet n, 

U^9A) = ^{(l-^i(r))( eVC o OS/i )+^i(r)f/(/i,0,0)(J)} 

Vd f e lfl (cos// — ihi(r) sin /x cos 6 1 ) \ to A\ 

72 V -e-^i(r) sin //sin^ J' ^ > 

*ufa,r,e,<l>) = ^{(l-/i 2 (r))( e _4 os/x )+^(r)C/(/.,0,0)(5 
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v u e 



ip 



^hiif) sin \x sin 6 



^[2 \e * M (cos/i + ih 2 (r) sin/icos#) / ' 



(2.5) 



n(/x,r,M) = ~^k(r). (2.6) 

Now the profile functions to be solved are /(r), hi(r), h,2(r) and /c(r). 
The lagrangian of the gauge-Higgs sector of the NMSSM is given by 

C = -\F; v F a ^ + (D^rf D^ d + D^ u + O^d^n - V ($ d , $ u , ra), (2.7) 
where the Higgs potential is given by 



Vb($ d , n) = m^^rf + m^t$ u + \n\ 2 - ( AA A n$ d $ u + ^A K n 3 + h.c. 



^2 + 9i /^t a . _ At a ^ 2 _l M - : 2 



+ |A| 2 |n| 2 + $t$ u ) + |A$A + /tn 2 | 2 , (2.8) 

The terms in the first line of (2.8) are the soft supersymmetry-breaking terms, those in 
the second line comes from the D-term potential, which are the same as those in the 
MSSM, and those in the last line are the F-terms, which are peculiar to the NMSSM. 
Here we adopt the same notation as in [8]. 

All the parameters in the potential with mass dimensions are to be determined once 
one specifies the breaking mechanism of the supersymmetry and the mass scale relevant to 
it. Instead of specifying supersymmetry-breaking mechanism, we express them in terms of 
the vacuum expectation values (VEVs) of the Higgs fields by requiring that the tree-level 
potential has vanishing first derivatives at the prescribed vacuum parametrized by the 
VEV of the doublet v , the ratio of the doublet VEVs tan (3 and the VEV of the singlet 
v n . As shown in [8] 5 , the derivatives of Vo with respect to the CP-even order parameters 
vanish at the vacuum, if the soft masses are given by 

m\ = ^ A -i^„^ n tan/3-im 2 w cos(2/5)-^L(t; 2 sin 2 /5 + ^), (2.9) 

(1 \ 1 I Al 2 

Rx - -n Vn j V n COt f3 + -77& 008(2/3) - LL. ( v * cos > (3 + vty , (2. 10) 

v 2 |A| 2 2 
m 2 N = (R x - Tlv n ) — sin (3 cos [3 + R K v n —v 2 -\k\ 2 vI, (2.11) 

V n 2 

while those with respect to the CP-odd order parameters vanish, if 

h = )fv n . (2.12) 



'Since we are working with g 1 = 0, we should set mz — mw in the results in [8]. 
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Here we defined the following parameters which are independent of phase convention, 
R x = -L R e(AA A e^)), I x = -^Im (\A x e^) , 
Rn = -^Re(^ K e 3 ^), 

K = Re (\K*e t(p - 2lp) ) , X = Im e i{p - 2ip) ) . (2.13) 
We use these relations to express the static energy functional which is 

E[$ d , $„, n] = J d 3 x [jF%F% + (D^ d y D& d + {D^D^ U + |^n| 2 + F ($ d , $ u , n)] , 

(2.14) 

in terms of the profile functions. Upon inserting the ansatz for the noncontractible loop, 
(2.1), (2.2), (2.4), (2.5) and (2.6), into (2.14), we obtain 



4lTV 



92 Jo 



d£ < 4 sin n 



I r/\2 i 2/ 2 (l f) 2 . 2 

(/ ) + 72 sm V 



cos 2 (3 sin 2 n 



eiK) 2 + 2{h\{\ -ff+a-hff 2 cos 2 /i 

-2^(1 - - /) cos 2 /i)" 



sin 2 /3 sin 2 /x 



+^ 2 sin 2 /3 cos 2 [(hj + h\- 2k(h 1 h 2 - 1) - 2) sin 2 /1 + (k — if 



em 2 + 2(h 2 (i - ff + (i - h 2 yf 2 cos 2 



-2/i 2 (l - /i 2 )/(l - /) cosV) 



+ 



N 2 e 



2 £-2 



C n N 2 i 



e 



sin /3 cos [2k{k- l)(/ii/i 2 - 1) sin 2 /i + (A; — l) 5 



+ |^[(/i 2 - l)cos 2 /3- (/i 2 - l)sin 2 /3] 2 sin 4 /i 
+ — sin 2 /3 cos 2 /3 sin 2 // cos 2 //(/ii — h 2 ) 2 



C £ 2 

\ — ^- sin 2 /3 cos 2 P sin 2 // 
2 - 1) [(hi - 

(A;2_l)2_^^(2A ; 3_3A ; 2 + l) 



(/ii/i 2 - I) 2 sin 2 fi (h 1 - h 2 ) 2 (2 + sin 2 //) 

3 



4 



(A; 2 - 1) [(hi - 1) cos 2 P+(hj- 1) sin 2 /3 

we,,.* ^ 2 ^ 2 e 2 



sin 2 /x 
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4nv 

92 



(c 2 sin 4 /i + ci sin 2 /i + c ) , 



(2.15) 



where £ = g2^o r and the primes on the profile functions denote the derivatives with respect 



to £. We introduced dimensionless parameters defined by 



A = 



m 



V r , 

N = — , 
K 



(Rx - nv n /2)v n 
rri^r sin (5 cos (5 rriyy 
R K N 



B = 



(2.16) 



Cn — 



gV 



c\ = 



|A| 2 


1 |2 
r _ l«l 




L-k - o 


gl ' 


#2 



and the coefficients in the last expressions are defined by 
8/ 2 (l - ff 



c 2 



/ « 

Jo 



- cos 2 /?/(l - - 2h x + fh) - sin 2 /3/(l - h 2 )(f - 2h 2 + fh 2 ) 



+ S ( ( ^ ~~ x) cos2 ^ ~~ (/ * 2 " x) sin2 /3 ) 2 ~ u sin2 ^ cos2 ^ hl ~~ h2) 



sin 2 /3 cos 2 /3 ((M 2 - l) 2 - hji x - h 2 f 



(2.17) 



Cl = 



jf di 4(f) 2 + ^ (e 2 (^) 2 + 2(/n - /) 2 ) + ^ W + 2(h 2 - /) 2 ) 



sin 2 /3 



+-£ 2 sin 2 /3 cos 2 (3 (h\ + h\ - 2k(h x h 2 - 1) - 2J 



H — ^— — — sin /3 cos /3 A;(A; — X){h\h 2 — 1) H 
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2C 

sm 2 /3cos 2 /3(/n-/i 2 ) : 



C^iV_ 2 e 2 
4 

^(A;') 2 + ^sin^cos^-l) 2 + 4 



■(A; 2 - 1) ([h\ - 1) cos 2 /? + (h\ - 1) sin 2 /?) 



Co 



/ da 2 



(2.18) 



sin/3cos/9(A; — 1)' 



(2.19) 



In the presence of the CP violation, the X-dependent term in the potential is proportional 
to sin n cos \x cos 0, which vanishes upon integration over 9 to yield the static energy of the 
configuration. This might suggest that an ansatz without spherical symmetry reduces the 
potential energy compared to the spherically symmetric one in the CP-violating case, but 
such an ansatz with less symmetry will increase the gradient energy. Hence we expect 
that the spherically symmetric ansatz yields the least energy configuration for each //, as 
long as the CP violation is not so large. 

The equations of motion for the sphaleron are derived by applying the variational 
method to the energy functional (2.15) with \x = n/2. After solving the equations, we 
check that the configuration is at least locally maximum at fi = n/2 along the noncon- 
tractible loop by calculating the coefficients c 2 and c±. That is, the coefficients calculated 
with the solution must satisfy 2c\ + c 2 > 0, for which E"(/i = n/2) < 0. 
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2.2 equations of motion 

The static energy functional at \i = n/2 is given by 
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iff + 



N 2 



i 



+|cos^[eW + 2(l-/) 2 ^ 



sin 2 pk 2 (h' 2 ) 2 + 2(1 -f) 2 h 2 2 



+- sin 2 (3 cos 2 (3 £ 2 (h\ + h 2 2 + k 2 - 2h x h 2 k - 1) 
8 



+ 



C n N 2 



sin cos /3£ 2 (2/c(A; - l)/^ - (A; 2 



+ ^[(/ i 2 -l)cos 2 /?-(/i 2 -l)sin 2 / 3 
+^£ 2 [sm 2 /3 cos 2 /3(/i 2 - l)(h 2 2 - 1) + iV 2 (A; 2 - l)(h\ cos 2 /3 + h\ sin 2 /3 - 1) 
+ ^e 2 (A: 2 - l) 2 - ^f?{k - l)(2k 2 - k - 1)1. 



(2.20) 



Applying the variational method to this functional, we obtain the equations of motion for 
the profile functions: 



no 



-/)(!- 2/)- 



(h\ cos 2 /3 + hi sin 2 /?) , 



2(1-/)% , A 



+ — sin 2 P(hi - h 2 k) + 



C n N 2 



tan (3 h 2 k(k — 1) 



(2.21) 



(2.22) 



+ h [[h\ - 1) cos 2 - (h 2 2 - 1) sin 2 p) + ^y ((hi - 1) sin 2 f3 + iV 2 (A; 2 - 1)) h u 



2(1-/)% , A 



e 



+ — cos 2 P(h 2 - hik) + 



C n N 2 



cot p h x k(k - 1) 



(2.23) 



h 



Jll [(^ _ i) C os 2 /3 - (h 2 2 - 1) sin 2 P) + y ((hi - 1) cos 2 P + N 2 (k 2 - 1)) h 2 

k"(0 + -^'(0 

A C 
= sin 2 p cos 2 P(k — h±h 2 ) + sin /3 cos P (h\h 2 (2k — 1) — k) 



C 



+ — (h\ cos 2 P + h 2 2 sin 2 P-l)k + C K N 2 (k 2 - l)k - Bk(k - 1). 



(2.24) 



In order to solve these equations, we must specify the boundary conditions for the profile 
functions. We relegate the derivation of asymptotic behaviors of the solution to the 



appendix and present only the result here. As shown in the appendix, the profile functions 
with finite energy functional satisfy the boundary conditions 



/(0)=/n(0)=/i 2 (0) = 0, fc'(0) = 0, 
/(oo) = /ii(oo) = h 2 (oo) = fc(oo) = 1. 



(2.25) 
(2.26) 



Note that only k(£) at £ = must satisfy the boundary condition of Neumann type, while 
the others satisfy those of Dirichlet type. 



In order to implement numerical study of the sphaleron solutions by solving the equations 
of motion (2.21) - (2.24) with the boundary conditions (2.25) and (2.26), we adopt the 
relaxation method which is suited for a boundary problem of ordinary differential equa- 
tions. The NMSSM has more parameters than the MSSM, but some of the parameters 
are constrained by the vacuum condition and the Higgs spectrum condition[8]. The model 
with the weak scale expectation value of the singlet scalar exhibits new features with re- 
spect to the Higgs spectrum and the phase transitions [7], which are absent in the MSSM. 
In particular, there are four types of phase transitions at finite temperatures. Although 
the sphaleron solution which is relevant to the decoupling of the sphaleron process at the 
first-order electroweak phase transition is a configuration which mediates the symmetric 
phase and the broken phase at the transition temperature, we here study sphaleron con- 
figurations with the tree-level potential for parameters satisfying the tree-level vacuum 
condition. 

In the case of the weak scale expectation value of the singlet scalar field, the allowed 
region in the parameter space is divided into two classes, one of which admits Higgs bosons 
lighter than 114GeV but with small couplings with the gauge bosons. We shall refer to 
this situation as the light-Higgs scenario. The other class of parameters yield heavy Higgs 
bosons satisfying the present bound on the neutral scalar. This situation is called the 
heavy-Higgs scenario, for which the Higgs spectrum and the phase transition are similar 
to those in the MSSM. 

We found sphaleron solutions for wide range of the allowed parameters, as long as the 
smallest eigenvalue of the mass-squared matrix of the CP-even Higgs bosons is positive, 
that is, the vacuum is at least a local minimum of the potential. This is obvious from the 
asymptotic behaviors of solutions, derived in the appendix. In the heavy-Higgs scenario, 
the energy of the sphalerons E satisfy 1.6 < E /(4nv / g 2 ) < 1.9 for wide range of tan/3 
and v n . The range of the sphaleron energy coincides with that in the MSSM. Further 
the energy becomes larger for larger mass of the lightest Higgs scalar. This behavior of 
the energy is expected from the results of the minimal standard model and those of the 
MSSM. As an illustration, we show the contour plots of the sphaleron energy and the mass 
of the lightest CP-even Higgs boson m Sl for tan f3 — 5, v n — 200GeV, m H ± = 600GeV and 



3 Numerical solutions 
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Figure 1: The sphaleron energy in the unit of 47ruo/(72 and the mass of the lightest CP- 
even Higgs boson as functions of A and k for tan/3 = 5, v n — 200GeV, m#± = 600GeV 
and A K = -lOOGeV. 

A K = — lOOGeV in Fig. 1. For this parameter set, both the heavy-Higgs and light-Higgs 
scenarios are realized for different regions in (A, n) -plane. Since we are working at the tree 
level, some of the parameter sets may be excluded by the Higgs mass bound. Our aim, 
however, is a qualitative study of the sphaleron solutions for a broad range of parameters. 
The shaded regions in Fig. 1, in which the smallest eigenvalue of the mass-squared matrix 
of the CP-even Higgs bosons is negative, indicates the parameters for which there is no 
sphaleron solution. We found that all the solutions satisfy the condition 2c2 + c\ > 0, 
which is necessary for them to be locally maxima along the noncontractible loop in the 
configuration space. These contour plots show that the sphaleron energy is larger for the 
larger mass of the lightest Higgs scalar in the heavy-Higgs regime (k < —0.5). 

In the light-Higgs scenario (—0.4 < k < 0), which is realized for small \k\, there 
is a minimum of the potential along the w n -axis. This minimum corresponds to the 
intermediate phase at finite temperature, at which the electroweak gauge symmetry is 
restored. Now we denote the potential in the symmetric phase (v = 0) as Uo(k) defined 
by U (k) = V (0,0,v n k/V2)/N 2 , that is, 

U (k) = °*f-k* - | k* + I (-A - C n - 2C X - AC K N 2 + 45) k 2 . (3.1) 

Depending on the parameters in the potential, Uo(k) has a nontrivial minimum at k 7^ 0. 
For smaller \k\, the value of k at the minimum becomes larger and the value of the potential 
at the minimum becomes smaller. Although the parameters which yield deeper potential 
at the minimum than the prescribed vacuum are excluded by the vacuum condition, some 
of the allowed parameters admit the new type of two-stage phase transition, which is 
referred to as type B phase transition in [7]. As seen from Fig. 1, the sphaleron energy 
is not so small in spite of the small mass of the lightest Higgs scalar, in contrast to the 
heavy-Higgs scenario. This small enhancement of the sphaleron energy is caused by the 
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Figure 2: Contour plots of the location of the minimum of Uo(k), k m m, and the value of 
k(0) of numerical solutions for the same parameter as those in Fig. 1. 



fc(0) 



-0.8 -0.6 -0.4 -0.2 



Figure 3: The values of k min and k(0) along the line of A = 0.85 in Fig. 2. 



discrepancy of the boundary value k(0) and the location of the minimum of the potential 
along the v n -sads. In Fig. 2, we show the values of k min at which Uo(k) is minimum and 
k(0) of numerical solutions in the (A, «)-plane The difference between k(0) and fc m i n is 
allowed because the boundary condition on fc(£) at £ = is of Neumann type, but not of 
Dirichlet type. The cross section of the contour plot at A = 0.85 is depicted in Fig. 3. The 
value of fc m in becomes larger than 100 for k ~ 0, while k(0) stays at a value near unity. 
The profile functions for k = —0.1, —0.5 and —1 along A = 0.85 are shown in Fig. 4. The 
values of the sphaleron energy of these profiles are within (1.84 — 1.86) x 4irvo/g2- For 
k = —0.1, fc m in ~ 2.99 so that the profile with fc(0) near this point lowers the value of 
the potential energy, while increases the gradient energy. Then the optimal profile which 
minimizes the total energy has k(0) ~ 1.08. In the heavy-Higgs scenario, where ~ 1 
for < £ < oo, the profile is similar to that in the MSSM. 
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Figure 4: The profile functions /(£) (solid curve), h\{£) (dashed curve), h<i{£) (dotted 
curve) and (dashed-dottted curve) for k = —0.1, —0.5 and —1 along A = 0.85 in 
Fig. 2. 

4 Discussions 

We have constructed the noncontractible loop in the configuration space of the NMSSM 
gauge-Higgs fields, from which we derived the equations of motion for the sphalerons and 
found a necessary condition for a solution to be a local maximum of the energy functional. 
We showed that the boundary condition of the singlet scalar at the origin is of Neumann 
type, in contrast to the other conditions of Dirichlet type. This difference in the boundary 
conditions in the Higgs profiles is due to the lack of gauge interaction of the singlet scalar 
field. 

We numerically solved the equations of motion and found solutions for wide range of 
allowed parameters of the model. The solutions in the heavy-Higgs scenario are similar 
to those in the MSSM, in the sense that stays almost constant and the energy of 
the solution takes the same value as that in the MSSM. The energy is larger for larger 
mass of the lightest Higgs scalar. This fact has been observed in the minimal standard 
model and in the MSSM. The sphaleron energies in the light-Higgs scenario are almost 
the same as those in the heavy-Higgs scenario, in spite of the small mass of the lightest 
Higgs scalar. This is because the difference between fc(0) and fcmin, at which the potential 
along v — is the minimum, enhances the potential energy. Although the Higgs potential 
of the NMSSM have negative contributions from the cubic terms, the sphaleron energy 
in the NMSSM is almost the same as in the MSSM, for wide range of parameters. The 
sphaleron energy is essential to determine the sphaleron decoupling condition after the 
electroweak phase transition. If we assume that the prefactor in the sphaleron transition 
rate in the NMSSM is the same as in the MSSM, the sphaleron decoupling condition 
is determined solely by the energy of the sphaleron 6 . Then the sphaleron decoupling 

6 The prefactor is composed of the zero-mode contributions and the factors coming from the Gaussian 
integrals of the positive modes. The zero-mode contributions are mainly determined by the global sym- 
metry of the solution, which is common to the MSSM and the NMSSM. The positive-mode contributions 
depend on the profile of the sphalerons. The effect of these prefactors on the transition rate is much 
smaller than the sphaleron energy on which the rate depends exponentially. 
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condition vc/Tc, where vc is the magnitude of the expectation value of the doublet Higgs 
field at the transition temperature T c , also applies to the NMSSM. 

All the numerical solutions we have found satisfy the condition 2c 2 + C\ > 0. As 
seen from the definitions of c 2 and ci, (2.17) and (2.18), a large negative contribution is 
expected to come from c\ when is far from unity. We could not find such solutions, but 
we cannot completely exclude the case of negative 2c 2 +Ci. Even if we encounter the case of 
negative 2c 2 + c\, it may imply that our ansatz is inadequate for such parameters. In that 
case, we will need to know the global structure of the energy functional to decide whether 
the sphaleron process is possible and whether it is suppressed after the electroweak phase 
transition. 

The sphaleron configurations obtained here are solutions with the tree-level potential. 
As shown in [5], the solutions with the finite-temperature effective potential differs in 
the energy from those with the tree-level potential by several percents in the MSSM. 
More precise decoupling condition of the sphaleron process in the broken phase will be 
obtained by solving the equations of motion with the effective potential with the radiative 
and finite-temperature corrections in our model. 
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A Asymptotic solutions 

We shall derive approximate solutions to the equations of motion for the sphaleron (2.21) 
- (2.24) in the asymptotic regions at £ ~ and £ ~ oo, respectively. For the energy 
functional (2.20) to have a finite value, the profile functions must behave as /(£) ~ £ 1 ^ 2+e/ , 
hi(0 ~ const. + ^~ l / 2+th i and ~ const. + £ _1 / 2 + e fc at £ ~ 0. Here e's are some 
positive constants. Similarly, all the profile functions must approach 1 as £ — > oo. More 
detailed behavior of them can be obtained by the approximate equations of motion in 
each asymptotic region. At £ ~ 0, (2.21), (2.22) and (2.23) are reduced to 



where we have used /(0) = 0. The solutions to the second and third equations are 
hi(0 £ an d hi(£) oc £~ 2 for % = 1 and 2, while the latter is excluded by the finiteness 
of the energy functional. Then (A.l) implies that /(£) oc £ 2 , if /(£) yield a finite energy. 
The asymptotic behavior of fc(£) is different from those of /ii(£) and /i 2 (£). This is because 



f"(0 * ^/(0 - \ (>i(0 cos 2 + hj(0 sin 2 (3) , 



MO, 



(A.l) 



(A.2) 



12 



the gauge interaction generates the dominant term proportional to 2/£ 2 in the right-hand 
sides of (A. 2), whose counterpart does not exist in (2.24). In the asymptotic region £ ~ 0, 
the equation (2.24) is approximated by 

k"(0 + |A/(0 = Aok® - Bk\i) + C K N 2 k 3 (0, (A.3) 

where we have put A = (A sin 2 (3 cos 2 (3/N 2 - 2C n sin (3 cos (3 - 2C X - 4C K iV 2 + 4B)/4. 
In order to solve (A.3), we set 

oo 

M0 = 7 + E«^ n+!/ , (a ^0,z/>0) (A.4) 

n=0 

where 7 is some finite constant. Inserting this expression into (A.3), we obtain, to the 
lowest order in £, 

v{y + l)a r _2 + (f + l)(f + 2)air _1 + • • • = (4, - B 7 + C K NW)l + 0(D- (A.5) 

If (A — £7 + C K N 2 ^ 2 )^ = 0, this equation holds only when v > 2. When (A — £7 + 
C K N 2 ^ 2 )^ 7^ 0, this is satisfied iiv — 2 and 

6a = (A - 5 7 + C K A 2 7 2 )7, ax = 0. (A.6) 

Although 7 and ao are not determined unambiguously at this order, it always holds that 
= 0) = 0. Therefore, we must impose the boundary conditions at £ = of Dirichlet 
type for /(£), ^i(£) and ^(O) an d that of Neumann type for 

At £ ~ 00, all the fields must approach unity for the energy integral to be finite. 
To study the asymptotic behavior, we put /(£) = 1 + A/(£), = 1 + Ahi(x) and 

fc(£) = 1 + Afc(£) in the equations of motion (2.21) - (2.24) and keep only the linear terms 
in the deviations from 1: 

A/"(0^ Q + |)a/~1a/(0, (A 

2 4 r 1 A?" 2 

A ^i(0 + 7i A/i 'i(0 - 7 sin 2 /3( A/11 - A/i 2 - Afc) + — — tan/3Ak 
£4 2 



+ 4 



i (A/i! cos 2 /3 - Ah 2 sin 2 /?) + C x (Ah 2 sin 2 /3 + N 2 Ak) 

A sin 2 /3 + cos 2 (3 . , f A + l „ \ • 2 „ a , 
^ ^ A/11 + f — + C A J sin 2 £A/i 2 

+ (-4^2 sin2 + ^ tan/? + C A ) N 2 Ak, (A 



2 4 f-o /V 2 

A/i''(f) + -=tA/i' 2 (0 ~ -cos 2 /3(A/i 2 - Ah! - Ak) + — — cot(3Ak 
^ 4 2 



-i (Ah! cos 2 /3 - A/i 2 sin 2 (3) + C A (A/i x cos 2 (3 + N 2 Ak} 
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A 



+ C x cos 2 /3A/ii + 



A cos 2 /3 + sin 2 (5 



Ah 2 



Ak"(0 + ^k'(0 



+ 
A 



A 



cos 2 /3 + — cot p + C x ) N' 2 Ak, 



AN 2 



(A.9) 



G 



■ sin 2 cos 2 P(Ak - Ah x - Ah 2 ) + ^ sin p cos P{Ah l + A/i 2 + Afc) 



+C A (A/i! cos 2 P + A/i 2 sin 2 P) + 2C K N 2 Ak - BAk 

._A sin 2 £ + ^ tan(3 + £^ cos 2 /?A/ll 



+ 



4AT2 



cos 2 /3 + ^ cot /3 + C x ^j sin 2 /3A/i 2 

/3 cos 2 P + sin p cos /3 + 2C K N 2 - B^j Ak. (A. 10) 



Now we introduce a 3 x 3 symmetric matrix P whose elements are 



Pn 

Pl3 
P23 
P33 



-(Asm 2 P + cos 2 p), pis 



+ Cx"j sin /3 cos ft, 



.^_ s in 2 /3 + ^tan/3 + C A ) .Yens i. 



P22 = -{Acos 2 P + sin 2 p), 



A 
AN 2 



G 



sin 2 P cos 2 P + — sin /3 cos p + 2C K N 2 - P. 



(A.H) 



The matrix P is related to the tree-level mass-squared matrix of the CP-even Higgs bosons 
AAg by P = A4 2 ^/(Am 2 v ). As long as we restrict the parameter sets to satisfy the vacuum 
condition, all the eigenvalues of P are positive definite. Then Eqs. (A. 8) - (A. 10) are 
written as 



/A/n(f) cos/T 
Ah 2 (Q sin/5 
V NAk(0 



A/n(0 cos /A 
P | Ah 2 {i)smp 
NAk(0 J 



(A.12) 



where we have defined A/ij(£) = £A/ij(£) and Ak(C,) = £A/c(£). If we denote the eigenvalue 
of P as A a , A/ij(£) and Afc(£) are expressed as linear combinations of e"^^"^ = e~ mSar , 
where ms a is the tree-level mass of the CP-even Higgs boson. Hence the asymptotic 
behaviors of the solutions are governed by the smallest mass of the CP-even Higgs boson. 
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